clear;
n=64;
k=0:0.01:n;
w=[1/3;1/2;2/3;1];
figure(1);
for i=1:4
y=1-2*w(i)*sin(k*pi/(2*n)).^2;
plot(k,y);
hold on
end
xlabel('k'),ylabel('y'),title('Eigenvalues of the iteration matrix R_ω'),grid on
saveas(1,'924Eigenvalues.png');